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Abstract 

Proton emission from deformed nuclei is described within the non-adiabatic weak coupling model 

O . 

which takes into account the coupling to 7-vibrations around the axially-symmetric shape. The 



coupled equations are derived within the Gamow state formalism. A new method, based on the 



^ ■ combination of the R-matrix theory and the oscillator expansion technique, is introduced that 

allows for a substantial increase of the number of coupled channels. As an example, we study the 
deformed proton emitter 141 Ho. 
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I. INTRODUCTION 



Theoretical models applied to the description of non-spherical proton emitters can be 
divided into two groups. The core-plus-particle models describe the radioactive parent 
nucleus in terms of a single proton interacting with a core (i.e., the daughter nucleus). 
Usually, the core is represented by some phenomenological collective model, e.g., the Bohr- 



Vfottelson (geometric) mode 



Depending on the structure of the daughter nucleus, rotational 



□ □ □ or vibrational 0, 0] couplings are assumed. The models belonging to this group 
employ the coupled-channel formalism of reaction theory which has been developed in the 
context of elastic or inelastic scattering. 

Models belonging to the second group employ the framework of the deformed shell model. 
In the simplest case, the proton resonance corresponds to a Nilsson state of a deformed mean 
field 0, Q, 0, |l0 ,^1 , Approaches belonging to this group can be generalized to include 
the BCS pairing lj 

We may refer to the first group of models as weak-coupling models or coupled-channel 
models. For the second group of models, we reserve the term resonance Nilsson-orbit (or 
adiabatic) models. The term "adiabatic" requires an explanation. It is very difficult to relate 
both groups of models to each other, because they operate on different approximation levels. 
In special situations, however, this relationship can be revealed. For instance, in the limit of 
the infinite moment of inertia of the axial weak-coupling model (which implies degenerate 
rotational bands and strong rotational coupling [14]), one recovers the resonance Nilsson- 
orbit model [16]. So one may say that in this case the adiabatic model is an approximation 
to the weak-coupling (non-adiabatic) picture. Generally, however, the relation between 
adiabatic and non-adiabatic descriptions is not simple. For example, the resonance Nilsson- 
orbit model with a triaxial potential jljj (i.e., nonzero 7 deformation) cannot be trivially 
related to a weak-coupling model extended to triaxial degrees of freedom jl^ . 

If the coupled-channel model with the rotational coupling is applied to the nucleus 141 Ho, 
the ground-state decay characteristics (half-life time and branching ratio) are poorly de- 
scribed There are several explanations possible. For example, it may be that the 
Coriolis mixing is too strong 3}. This can be partly cured if pairing is introduced 
Another possibility, explored in this work, is the coupling to triaxial vibrations. Indeed, 
in particle-plus-rotor calculations, the best description of the experimentally observed band 



2 



structures of 141 Ho can be explained if 7 deformation is considered |l7J]. In addition, in 



the neighboring nuclei, such as Sm and Gd, there are low-lying 1\ and 3 + levels 
which have been interpreted Il9j as members of a 7-vibrational band. There are also other 

nn 

indications that in this mass region the coupling to triaxial modes can play a role [20. 12 1|. 
The possibility that triaxiality influences the decay of 141 Ho was investigated in our earlier 
work Q| and also in the recent Refs. 15, based on an adiabatic model assuming a 



triaxially deformed mean field. 

In this work, we present non-adiabatic calculations in which t re excitations of the daugh- 
ter nucleus are properly taken into account. Unlike in Ref. 22], we do not assume a per- 
manent 7 deformation of the core, but rather we consider 7 vibrations around the axially- 
symmetric deformed shape. 

The ground-state rotational band of 140 Dy has recently been observed jQ, .24] . In addi- 



tion, in our work we assume that 140 Dy has the K=2 7-vibrational band. This structure can 
be coupled to the ground-state band if the proton-daughter interaction in the body-fixed 
system deviates from the axial symmetry. The experimentally observed rotational band of 
the parent nucleus is assumed to be a K 7T =7/2~ band built upon the [523]q =7 / 2 Nils- 
son level. In the strong-coupling picture the presence of the 7 band in 140 Dy implies the 
existence of two additional rotational bands in 141 Ho with K=Q ± 2, i.e., = 3/2" and 
K* = 11/2- 

In the weak-coupling model, proton emission is described by means of a coupled set of 
differential equations which are solved assuming appropriate boundary conditions. The most 
obvious way to describe the proton emission is to assume outgoing boundary conditions. This 
immediately leads to the notion of the Gamow or resonant states, the generalized eigenstates 
of the time-independent Schrodinger equation, which are regular at the origin and satisfy 
purely outgoing boundary conditions. Together with non-resonant scattering states, Gamow 



states form a complete set, the so-called Berggren ensemble |25j, which can be used in a 

n nnn 

variety of applications 26], including the recently developed Gamow shell model 27, 28, 29] . 



Unfortunately, the number of coupled equations rapidly increases with the number of 
excited states of the daughter nucleus taken into account. In addition, the solution of the 
eigenvalue problem of a very large set of coupled equations becomes numerically unstable 
at some point. This is especially true if one keeps in mind that there is a twenty-order-of- 
magnitude difference between the real and imaginary part of the energy of the Gamow-state 
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which describes the proton decay of 141 Ho. A possible way out is to consider the R-matrix 
theory. However, even in this case, one has to deal with large sets of coupled differential 
equations. 

In order to avoid the difficulty of solving large sets of coupled differential equations, one 
may use the Rayleigh-Ritz variational principle and apply the basis expansion method. In 
this paper, the spherical harmonic oscillator wave functions are used as basis functions. It 
was recognized a long time ago that by using the basis expansion method the positions of 
narrow resonances can be determined. In particular, the signature of a narrow resonance is 



that the specific 



positive energy 



solution is locally stable with respect to the change of the 



size of the basis 3(], 31, 32, 33, 34, 35, 36,|s7|. Several proposals exist in the literature on 
how to determine the width of the resonance in this method. They are called L 2 stabilization 
methods |34j . (The name comes from the fact that only square integrable functions are used 
in the expansion.) In this paper we will introduce a new method which is a combination of 
the oscillator expansion method and the R-matrix formalism. This method is very simple 
and proves to be accurate enough for very narrow proton resonances. 

The paper is organized as follows. We begin in Sec. |TI] with an overview of the weak- 
coupling model applied to the case of rotational motion and 7 vibrations. Sect ion UTTl reviews 
different methods to calculate the position and width of a resonance state: the theory of 
Gamow-states, the standard R-matrix formalism, and the new method which combines the 
oscillator expansion method with the R-matrix formalism. Finally, Sec. IIVI contains results 
of numerical calculations. We check the accuracy of the new method and demonstrate how 
the position of excited states in the daughter nucleus can influence predictions of the weak- 
coupling model. We also present results for the proton emission in 141 Ho. Finally, Sec. IS 
contains the conclusions of this work. 

II. WEAK-COUPLING MODEL 

The proton-emitting parent nucleus is described here in terms of a single proton coupled 
to a deformed core. The model Hamiltonian can be written as 

#rot = H d -^-A r +V dc{ (r,u), (1) 
2m 
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where is the (collective) Hamiltonian of the daughter nucleus, the second term represents 
the relative proton-daughter kinetic energy, and Vdef is the proton-core interaction, which 
depends on the position of the proton r and the orientation u of the core. 



A. The proton-daughter interaction 

It is straightforward to define Vd e f in the body-fixed frame, in which one can define 
the deformed mean field. By expanding the nuclear radius in multipoles and assuming 
quadrupole deformations only, one obtains 

R(9', 0') = R C(a , a 2 ) [1 + a Y 2>0 (9') + a 2 (Y 2j2 (9', 0') + F 2) _ 2 (0', cj>'))\ (2) 

where C(ao, a 2 ) is the volume conservation factor. The intrinsic deformed field is defined 
using a Saxon- Woods form factor 

V dei (r, m = -— V ° (3) 

1 + exp [(r — R(9, 4> ))/a\ 

Expanding to the first order in a 2 , one obtains 

Vdef (r, 9'4>') = V x {r, 9') + a 2 V 2 (r, 9') [Y 2)2 (e', 4>') + Y 2 ^ 2 (9', </>')] . (4) 

The form factor Vi(r,9') is the same as Q except that a 2 is put equal to zero. The form 
factor of the second term is given by 

V 2 {r,9) = (5) 

a 1 + e 

where, again, a 2 =0 in R(9', <f/). The deformed form factors V\{r, 9') and V 2 (r, 9') still depend 
on ao- After perfoming multipole decomposition of V\ and V 2 , one obtains the intrinsic 
potential: 

V def (r, 9'<f/) = V$(r, 9') + a 2 v£l(r, 9><f>>) 
= Ea V^(r)Y{ fi (9') + a 2 £ A V A (2) (r) [Y^ff, 0') + 0')] • (6) 

For explicit expressions for V^\r) and V^\r) see, e.g., Ref. It can be shown that in 
the laboratory system the daughter-proton interaction is given by 

V def (r,u,) = V^(r,w) + OaVjJfcw) 
= Ea, V A (1) (r)^ o r A , M (f) + a 2 Ea, Vf(r) (D^ + D^ 2 ) Y x ,,(r). (7) 
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In addition to the nuclear potential, there is also a long-range Coulomb interaction between 
the deformed core and the proton. The deformed Coulomb form factors, Vq and Vq \ are 
discussed in Appendix A. 



B. The coupled channel equations 



The states of the daughter nucleus are eigenvectors of In this work, we adopt the 
rotational- vibrationa 
the standard ansatz 



collective model. The wave functions of the core, (pi^x, are given by 
3: 



h,K = J 1Q J(s^+ 1) + ( -1 J Wa 2 )|g.s.>, (8) 

where XKn 2 { a 2) is a 7- vibrational wave function. The wave function of the parent nucleus 
can be written in the weak-coupling form 

&m = J2 u j^l^ Mmi , (9) 

where the channel function is given by 



r 

IKlj 



and 



QjMiKij = ^(jm^\jM)y m (j) I ^ (10) 



y m = Y,(lrnUjn)t l Y lm (r)xi/2(s) (11) 



ms 

arises from the coupling of the proton spin with the orbital angular momentum. In our 
earlier weak-coupling calculations there was no summation over K in Eq. JOJ); only the 
K — term was considered. Due to the non-axial symmetric form of the proton-daughter 
interaction (|7|). the ground state K = and the 7- vibrational K = 2 band both contribute. 
The radial functions uj Kl j(r) are solutions of the set of coupled-channel equations: 

£ (-£* + ^) u J IKlj + Exrvr MHj, ri'f, J)B x {II'K)vPu J rKVjl + (12) 
Exi>K>i>j> Mty, J)Cx(IKI'K', a 2 )V^u J rKIVjl = (E - E IK )u J IKlj , 

where Ejk is the energy of the daughter state described by the wave function (jHJ). The 
r-independent coupling coefficients can be written in terms of the reduced nuclear matrix 
elements 

B X (II'K) = (MID^pk) (13) 



and 



C x (IKI'K',a 2 ) = (<MMD A 2 + D A _ 2 )||0^>- 



(14) 



The explicit expressions for the geometric coefficients A\(Ilj, I'l'f , J) are given, e.g., in 



Ref. 



. The nuclear structure model of the daughter nucleus enters the formalism through 



the reduced matrix elements B\ and C\ 



38, 



III. CALCULATION OF RESONANCE PARAMETERS 

The coupled differential equations (|12j) can be turned into an eigenvalue problem by 
specifying boundary conditions. It is always assumed that the solutions are regular at the 
origin, i.e., u c (0) = 0. (From now on, the channel indexes IKlj are abbreviated by the 
symbol c.) 



A. Gamow states 



To be a Gamow state, the radial wave function must asymptotically behave as an outgoing 
Coulomb wave: 

u c {r) Oi(r],rk p ) 

= Gi(r),rk p ) +iFi(rj,rk p ), (15) 

where k\ = j?r(£ p — Ejk) and rjk c = j^Ze 2 . Such boundary conditions are only satisfied for 
a discrete set of complex wave numbers k c which define the generalized eigenvalues E = £ p 
of Eq. (|12|). These eigenvalues correspond to the poles of the scattering matrix 

Q,Q- The 

corresponding solutions are either bound states with negative real energies S p = < and 
pure imaginary wave numbers k p = (j p > 0), or resonance states, £ p = E TCS — zEff% with 
nonzero imaginary parts r rcs ^ 0, and k p = k p — ij p . 

The asymptotic behavior of the radial wave functions are determined by k p . For Gamow 
states these functions show oscillating behavior at large values of r so one must define a new 

n 

normalization scheme. Berggren proposed [25| a generalized scalar product and introduced 
a regularization procedure (IZeg). With this generalization the norm is 

POO 

J2 ne 9 [u c {r)fdr = l. (16) 



Once the resonance energy and radial wave function have been determined, there are 
different methods to calculate the width of the state. The simplest method is to take twice 
the imaginary part of the energy of the resonance. However, for narrow resonances the 
accurate numerical calculation of Zm[£ p ] is difficult. Therefore, other methods are often 
used. One possibility is to calculate the partial width for each channel from the so-called 
current expression [40] 

fi 2 u'*{r)u c {r) - u' c (r)u*(r) 
C{) 2„ £ c J r k'(-')| 2 dr' ' li7J 

where the sum of the partial widths 

r res = ^r c (r) (is) 

c 

gives the total decay width. Although values of T c (r) depend on r in the region where the 
coupling potential terms are not negligible, the total width (|T8|) is independent of r, which 
reflects flux conservation. 

In practice, the Gamow boundary condition given by Eq. (fTH|) can be implemented in the 
form 

<( r as) = k 0[{r] ) r as k p ) 

^c(^as 

where r as is the channel radius (the off-diagonal couplings are negligible for r > r as ). Using 
Eq. (|19J) . the partial decay width can be written at the point r as as 

r c (r as ) : 



h 2 


Mc^as)! 2 


2/z lOjfofcpras) 


l 2 E c j; aS M^)l 2 dr' 



X 



[k*0'*(r), r as fc p )0;(r7, r^kp) - k p O[(r], rjz p )0*{r], r as /c p )] . (20) 



If one neglects the imaginary part of k p , the square bracket in Eq. (J20|) becomes — 2i and 
the expression for the partial decay width can be written in a simple form: 



r^as) Hi ^ >,,^ ffas I I A 19 J /' (^1) 









Oi(rj, fc p r as )| 


2 V- fas 
Z^c' Jo 


u c '(r')\ 2 dr' 



Equation (J2Uj) and its approximate form (|21|) are strictly valid only at the point r as where 
the boundary condition is given. We emphasize at this point that if the coupled equations 
are solved with the Gamow boundary condition, then the total width can be calculated at 
any value of r using exact relations (fTTjl and IjIBj) . 
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B. R-matrix method 

For completeness, we summarize those important aspects of the R-matrix theory [41] 
which are relevant to our work. In the R-matrix theory one also deals with a set of radial 
functions g c (r). These functions are regular at the origin and satisfy the same coupled 
equations (fT2"j) as the Gamow states but with the following boundary conditions 

M = B e , (22) 
9c{a) 

where the parameters B c are arbitrary real numbers. It is assumed that the short-range 
diagonal and off-diagonal proton-core interactions can be neglected beyond the channel 
radius a. Consequently, a has the same meaning as the parameter r as of the Gamow theory. 
It is worth noting, however, that a is always real, while r as can be complex. 

The boundary condition (}2*2*j) defines the complete set of functions inside the channel 
surface. The real eigenvalues of the coupled-channel equations are denoted by E x and 
the corresponding eigenfunctions by (r). They are normalized to one inside the channel 
surface, 

£ f \9c(r)\ 2 dr = 1, (23) 

c J 

and define the so-called reduced width amplitudes 

/ h 2 \ 1/2 

The resulting R-matrix has a simple form 

^) = E|fr§ (25) 



but it is related to the physically important scattering S-matrix in a complicated way (41 1. 
Let us emphasize that the calculated S-matrix is independent from both the boundary 
condition parameters B c and from the channel radius a only if all the R-matrix states are 
taken into account in Eq. (|25|). 

Assuming that in a given energy region only one term dominates in the R-matrix and 
making further approximations (see p. 322 of Ref. jl^), Lane and Thomas showed that the 
S-matrix can be written in the form 

iY Xc (E) 1/2 T Xd (E) l/2 



S CC ,(E) « S° CC ,(E) 



0~c.c' + 



E x + A x {E)-E-ir x {E) 



(26) 



where the partial R-matrix widths 

T Xc (E)=2P lc (Ehl (27) 

give the total width 

T x (E) = J2^Xc(E). (28) 

c 

In Eq. (J2H1), function A X (E) is given by 

A A (£) = ]TA Ac (£), (29) 

c 

where 



A Xc (E) = -(S lc (E)-B c ) 1 l. (30) 



The penetration Pi c (E) and shi 



functions (see p. 270 of Ref. |4lj). 



t Si c (E) functions are related to the Coulomb Fi c and Gi 



Within approximation (J26)) . the complex-energy resonance poles of the S-matrix, 
— |r^ s , satisfy the equation 

■>R\ cR ! r (cR 



E x + A A (5«) - f* - -T X (E«) = 0. (31) 

Here, we used the upper index R in order to distinguish this R-matrix approximation for the 
resonance energy from the energy of the corresponding Gamow state. In order to simplify 
the solution of the non-linear equation (|3ip. one often introduces further approximations 
and assumptions for the calculation of the functions A X (E) and T X (E). 



42|, the function A Xc (E) is expanded around the R-matrix 



In the method of Thomas 
eigenvalue E x 

A Xc (E) w —(Si c (E) - B c ) - S lc (E x )(E - E x ), (32) 

where the dot denotes energy derivative. Furthermore, the ^-dependence of T Xc (E) is neg- 
lected and T Xc (E) is replaced by the corresponding value at E x . Under these assumptions 
one obtains 

e r = E x + ^ c (B c -S lc (E x )) 7 l 
l + EMExhl 

and 

r R = Tx{Ex) . (34) 
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In order to simplify ()33|) we may require that the chosen boundary condition parameters 
satisfy the condition: 

B c = S lc (E x ). (35) 

If the Si c (E\) terms are negligible, then the resonance energy corresponds to the R- matrix 
eigenvalue 

= E x (36) 
and the width can be calculated with the well-known expression 

T^ s = J2^Pi c (Exhl c . (37) 

c 

Two variants of Thomas's procedure can be found in a later paper of Lane and Thomas 



4l| where they give different expressions for E^ and T^ s . 



C. R- matrix method using oscillator expansion 

In this section we propose a simple method, based on the R-matrix formalism, to estimate 
the parameters of a resonance. The advantage of this method is that it avoids solving a large 
set of coupled differential equations. The method is based on the expansion of the radial 
functions ujKij{r) in the single-particle basis 4>ni°( r ) °f the spherical harmonic oscillator. In 
this basis, the total wave function (jUJ) can be written in the form: 

* JM = E E Ci K J-^^ JMI K lr (38) 

IKlj n 

The coefficients Cj Kn[ ^ can be obtained from the matrix eigenvalue equation: 

/aHO\H 2 ( d 2 , KH-lA \iHO\ryJ I r?HO r? \r<J 

l^n'\ ( Pnl I 2^ ydr 2 ^ r 2 J Wn'l I^IKn'lj ~ ~ ^IKj^IKnlj 

+ Eaw MW, Mf, ^)5A(//^)(0nn (1) l«)C/'^Y (39) 
+ Expkw MHO, I'l'f, J)C X (IKI>K>, a 2 )(^°\vi 2) \<P^)Cl Kfnll , r = 0. 

In the following, the corresponding real eigenvalues are denoted as E x ° . 

In the R-matrix theory, the coupled equations (fT2j) are solved with imposed boundary 
conditions However, as discussed in the following, this procedure can be reversed. In 

the first step, we solve the algebraic eigenvalue problem (|3TI|) for the coefficients Cj Knl y The 
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resulting radial functions g c {r) define the boundary condition function at the point r: 



\ n / n 



B c (r) = B IKlj {r) = r ( £ C{ Knlj <j>^{r) ) / £ C/ Kn ^°(r). (40) 



Having determined the boundary condition parameter at each r, the R-matrix formalism 
can now be applied. In particular, after replacing E\ with in expressions and (pMj) . 
they can be used to compute the position and the width of a resonance at each value of r: 



£f° + E c (Se(r) - S^E* °)) 7Ac 



r 



,2 



and 



resU l + E c ^(^°)TA C (r) 2 1 J 

r HO M = r A(£f°) C42I 

resU l + E c ^W°)7A C (r) 2 ' 1 J 

where the r-dependent reduced width amplitudes (|24|) are given by 

7Ac(r)=(^) 1/2 $:C/^ (r). (43) 

This algorithm is further referred to as the R-matrix method based on harmonic oscillator 
expansion (RMHO). In RMHO, the energy and width of the resonance explicitly depend on 
r. However, for sufficiently large values of r, this dependence is expected to be extremely 
weak. It is to be noted that since expression (J3*T|) is derived under specific assumption (|3*3jl . 
it is not valid in the RMHO method. 

The derived boundary condition parameters (j4*Uj) do not depend on the actual normaliza- 
tion used. However, this is no longer true for the reduced width amplitudes (j24|) . In order 
to apply the R-matrix method at each a=r, the radial functions 



'■max 



9c(r) = gjKiAr) = £ Cj Knl ^{r) (44) 



n=0 



have to be renormalized to one inside the channel surface according to Eq. (J23)) . 



IV. RESULTS 

The numerical tests have been carried out for the deformed proton emitter 141 Ho, viewed 
as a proton-plus-core system, with the daughter nucleus 140 Dy being the collective core. We 
employed the same successful parameterization of the Woods-Saxon (WS) optical potential 
as in earlier Ref. Q]. 
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A. Resonance Width in RMHO 

Let us first assume that the core is axially deformed (ei2 = 0). In the calculations, all the 
states in the g.s. rotational band in the daughter nucleus up to 1=12 were considered. In 
our weak-coupling calculations, the experimental excitation energies of 140 Dy were used for 
states with J<10, and the energies of the remaining states were obtained by the variable- 
moment-of-inertia (VMI) fit to the data. That is, for the g.s. band we took the values: 0.203, 
0.567, 1.044, 1.597, 2.218, and 2.894 MeV. The deformation parameter ao was set to the 
value of 0.244, which is consistent with earlier investigations jl 71. |23|. The WS strength was 
adjusted to reproduce the experimental position of the J 7r =7/2~ resonance at 1.19 MeV. 
The number of coupled channels in this variant is 46. This number is sufficiently small 
to carry out the reliable calculation of the Gamow-state energy eigenvalue. The resulting 
resonance width is 0.208 x 10 -19 MeV. We accept this number as the exact, or reference, 
value. 

The harmonic oscillator basis is characterized by a single parameter, the oscillator length 
b. The upper part of Fig. Q shows the resonance width (J42|) calculated in RMHO as a 
function of r. For each partial wave, M=n max +l=12 harmonic oscillator functions were 
used in the expansion (J4*4^ and the value of b was varied. As expected, a clear plateau 
appears at large values of r. The extent of the plateau depends on the size of b: the greater 
oscillator length (i.e., the r.m.s. oscillator radius), the greater the extent of the plateau. The 
reason for the rapid decrease of the width function r^f (r) at very large values of r lies in 
the fact that the radial channel function is approximated by a linear combination of a finite 
number of oscillator functions, each having the Gaussian asymptotic behavior. Therefore, 
by increasing the number of states in the basis, the extent of the plateau is expected to 
increase. This is illustrated in Fig. ^ (lower portion) which shows RMHO results obtained 
at a fixed value of 6=1.8 fm for several values of M. It is seen that for M=24 (n max =23) 
the width function becomes independent of r in a very wide interval of r. In the interval 
between r=9 and 12 fm the RMHO width exhibits tiny oscillations (practically invisible 
in Fig. [IJ. Therefore, to obtain a well-defined value, we divide this interval equidistantly 
with a step size of 0.1 fm and calculate the average F^P = ^- (?"«), which will be 

considered as the RMHO width in the following. 

In order to assess the quality of the RMHO method, Fig. |2] shows the relative errors 
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FIG. 1: The width of the J*=7/2~ resonance in 141 Ho at 1.19 MeV calculated in RMHO as a 
function of r. Top: dependence on the oscillator length parameter b (the number of basis states is 
M=12). Bottom: dependence on M (6=1.8 fm). 

of the real and imaginary part of the energy of the resonance as a function of the WS 
potential depth V . The reference values were obtained by the Gamow-state coupled-channel 
procedure. In the considered region of Vo, the resonance width changes by four orders of 
magnitude; however, the relative error of RMHO is less than 1.7 percent. The accuracy of 
RMHO for the real part of the energy is much better: the relative error is always smaller 
than 0.0025 percent. The results presented in Figs. ^ and El convincingly demonstrate that 
the RMHO formalism can be safely used to calculate isolated narrow proton resonances. 
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FIG. 2: The relative error of the real (solid line) and imaginary (dashed line) energy of the J w =7/2~ 
resonance in 141 Ho calculated in RMHO (M=24) as a function of the WS potential depth Vq. The 
reference (exact) values are taken from the Gamow states calculation. In the considered range of 
Vq, the resonance width changes by four orders of magnitude. Note that the solid curve has been 
multiplied by a factor of 100. 

B. Proton decay of 141 Ho 

In this section we investigate the influence of 7 vibrations on the process of proton 
emission from 141 Ho. All results presented in this section are obtained with the RMHO 
method using M=20 oscillator functions for each partial wave. The oscillator length was 
assumed to be 6=1.8 fm. Using the results of the VMI fit for the g.s. band, the assumed 
energies of the members of the 7 band are: 0.750, 0.934, 1.144, 1,378, 1.633, 1.907, 2.198, 
2.504, 2.825 3.159, and 3.507 MeV for I = 2,3,4, ...12. The chosen position of the 2^ 
band head of the 7- vibrational band, 750 keV, was taken according to the systematic trends 
around iV=74. 

When one includes the K=2 7-vibrational band in addition to the g.s. band, the number 
of coupled channels increases from 46 to 130 (assuming that the maximum spin is J max =12 
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TABLE I: $7 decomposition of the J n =7/2 states in 141 Ho in the energy region (-3 MeV, 1 MeV) 
calculated in the non-adiabatic approach. The axial deformation (ao=0.244, a2=0) is assumed. 



E r (MeV) 


n=l/2 n=3/2 17=5/2 n=7/2 


-2.255 
-1.066 
-0.103 
1.190 


0.828 0.163 0.009 0.000 
0.168 0.694 0.135 0.003 
0.006 0.144 0.808 0.042 
0.000 0.001 0.045 0.954 



in both bands). For that reason, we decided to carry out the RMHO calculations instead of 
the weak-coupling Gamow analysis. We have checked, however, that for / max = 10, where the 
coupled-channel calculations with the K=2 band can be done, the RMHO results coincide 
with those of the coupled- channel method. 



1. Structure of J 77 =7/2" states 

The simplest model of the g.s. decay of 141 Ho is based on the adiabatic resonant Nilsson- 
orbit picture of Sec. II. B of Ref. Q]. Here, the valence proton occupies a Q 7T =7/2~ Gamow 
state in an axially deformed mean field. Let us consider this scenario first. In our WS model, 
there is only one Q=7/2~ state in the energy region (—3 MeV, 1 MeV) and deformation 
(a ~0.244, a 2 =0). The calculated energy of this [523] 7/2 state is 0.426 MeV. There are 
three more negative parity Nilsson states originating from the /in/2 proton intruder shell, 
with energies -2.678 MeV ([550]l/2), -2.141 MeV ([541]3/2), and -1.105 MeV ([532]5/2). 
If we now apply the weak-coupling model (i.e., we assume that the daughter nucleus has a 
g.s. rotational band with the finite moment of inertia), we calculate one J n =l/2~ state, two 
3/2~ states, three 5/2~ states, and four 7/2 states in the considered energy region. Of 
those four 7/2~ states, only one can be associated with the 7/2~ band-head from which the 
proton emission takes place. The remaining three are rotational excitations associated with 
the K paiT =Q=l/2, 3/2, and 5/2 bands built upon the deformed Nilsson levels mentioned 
above. Table |J displays the structure of the J n =7/2~ states calculated in the non-adiabatic 
approach. The Q decomposition of the states Q clearly identifies the Nilsson orbit upon 
which the rotational g.s. band of the parent nucleus is built. 

The situation becomes more complex if, in addition to the g.s. band, one also considers 
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TABLE II: Same as in Table [I] but in the presence of small triaxial coupling {a 2 = 0.05). 







K 


=0 




K 


=2 


(MeV) 


0=1/2 


0=3/2 


0=5/2 


0=7/2 


0=1/2 


0=3/2 


-2.616 


0.691 


0.232 


0.029 


0.007 


0.034 


0.007 


-1.184 


0.015 


0.323 


0.004 


0.014 


0.429 


0.215 


0.122 


0.028 


0.078 


0.529 


0.214 


0.134 


0.017 


1.153 


0.015 


0.003 


0.018 


0.902 


0.002 


0.060 



the K=2 rotational band in the daughter nucleus (i.e., if one takes the non-zero triaxial 
coupling a 2 ). The coupling to 7 vibrations immediately results in an increase of the numbers 
of predicted bands. Indeed, since the 7 band can be built upon each K pauT =Q structure, one 
obtains twelve bands with quantum numbers K par ,, K pw , + 2, and \K paT , — 2| in the energy- 
interval considered. Among those twelve bands, only two have a J n =7/2~ band head. One 
is the previously discussed [523] 7/2 band while the other corresponds to a 7-phonon built 
upon the [541]3/2 Nilsson orbital. Table ILT1 displays the Q decomposition of the four states 
of Table U in the presence of a small 7 coupling (a2=0.05). It is seen that the single-proton 
band head is clearly identified. 



2. Proton emission from the ground state of Ho 

Earlier investigations have demonstrated that in the weak coupling model there is 

a sensitivity of the resonance's parameters to the number of states in the rotational bands of 
the daughter nucleus taken into account. Figure El shows calculated energies of the J 7r =7/2~ 
states in 141 Ho as a function of the coupling constant a 2 (a Q = 0.244) for several values of 

-^max- 

Figure El shows that for bound states the convergence is already very satisfactory for 
-fmax=10; however, this is not true for the 7/2~ band head. For instance, at a2=0.1 one 
obtains for the energy of lowest bound state: -3.254, -3.264, -3.265, -3.265, and -3.265 MeV 
for / max =8, 10, 12, 14, and 16, respectively. In contrast, for the 7/2" g.s. (marked by thicker 
lines in Fig. EI) the analogous numbers are: 1.297, 1.171, 1.062, 1.062, and 1.062 MeV. That 
is, in this case, going from J max =10 to J max =12 the energy changes by as much as 109 keV. 
This variation is significant since the width of the resonance is extremely sensitive to its 
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FIG. 3: The position of the bound and resonance J 7T =7/2~ states in 141 Ho calculated in the weak 
coupling model as a function of the triaxial coupling constant 02 for I max =10 and 12 (I max is the 
maximum value of angular momentum considered in the g.s. band and in the gamma band of the 
daughter nucleus 140 Dy). The results for J max =12 are fully converged, i.e., a further increase of 
the number of states does not change results in the scale of this figure. The axial results without 
coupling the 7 band (cf. Table. [IJ are marked by the dots. The 7/2~ ground-state of 141 Ho is 
marked by thick lines. 

energy. In our previous paper |16j, we made the pilot studies of the coupling to the 7 band 
on proton emission in 141 Ho. Unfortunately, in this early analysis based on the coupled- 
channel method, we took J max =10; hence the conclusions of this paper have to be revised 
(see below). 

The width of the 7/2" band head was computed using the RMHO method assuming 
I max =12. At each value of a 2 we have adjusted the potential depth so as to get the position 
of the resonance at 1.19 MeV. The calculated half-life of the resonance and the branching 
ratio for the decay to the 2f state in 140 Dy are displayed in Figure 0] It is seen that when 
increasing the coupling to the 7-band, both the lifetime and the 2± branching ratio increase, 
and the agreement with experiment gets worse. 
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FIG. 4: Half- life (top) and branching ratio for the decay to the 2f state in 140 Dy (bottom) for the 
7/2~ g.s. of 141 Ho as a function of the triaxial coupling a 2 . Experimental data (7^=4 ms and 
branching ratio 0.7%) are taken from Refs. 43 1. The spectroscopic factor (BCS occupation 
coefficient) was assumed to be w 2 =0.84 

In order to understand the behavior shown in Fig. we analyzed the components of the 
wave function. Figure El shows the weights of various partial waves (IKlj) in the 7/2~ g.s. 
of 141 Ho, 

POO 

\CiKij\ 2 = / uj Klj (r)dr, (45) 
Jo 

as functions of a 2 . According to our calculations, the amplitudes associated with the coupling 
to the 0^ g.s. and 2f state in 140 Dy are fairly small; most of the strength lies in higher- 
lying states including the channels that are energetically closed for proton emission. The 
(0+, / 7 / 2 ) amplitude, solely determining the 7/2~— >0f decay, gradually decreases with a 2 . 
Interestingly, while the total strength increases with a 2 as expected (the hn/2 and fj/2 
shells are strongly coupled by triaxial field) , most of this strength is pushed up to higher-lying 
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FIG. 5: Weights (gSJ) of the (J=2i, 1=5 j=ll/2), (Oi, 3, 7/2), (2i, 5, 9/2), and (2i, 1, 3/2) partial 
waves of the 7/2~ g.s. of 141 Ho as a function of 02- 

states. 

Figure El displays partial widths r c (|T7|) corresponding to various channels of decay to 
the Of and 2+ states in 140 Dy. The gradual decrease of the (Oi, /r/2) partial width (hence 
the increase of the half-life of 141 Ho) with a<i can be explained in terms of the (Oi, iV/2) 
amplitude in Fig. |3J The 2 + branching ratio is almost completely determined by the (2 1; 
jfV/a) partial width; the second-order contribution from the P3/2 wave is much smaller (cf. 
inset in Fig. EJ) . 

V. CONCLUSIONS 

This work contains the first application of the triaxial non-adiabatic weak coupling ap- 
proach to the description of proton-emitting nuclei. The resulting coupled-channel equations 
take into account the coupling to the K=2 band representing collective 7 vibrations. 

The inclusion of the 7 band into the weak-coupling formalism increases the number of 
the coupled channel equations significantly. This makes it very difficult to solve accurately 
the multitude of coupled differential equations with Gamow boundary condition. In order to 
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le-26 : 



1=0 1=3 j=7/2 
1=2 1=1 j=3/2 
1=2 1=3 j=5/2 
1=2 1=3 j=7/2 
1=2 I=5j=9/2 
1=2 1=5 j=l 1/2 

I 

0.04 



0.08 



0.12 




0.16 



FIG. 6: Partial widths r c (|17|) corresponding to various Of and 2f decay channels as functions of 
a 2 . 

overcome this difficulty, we developed a new formalism, dubbed RMHO, which incorporates 
the variational oscillator expansion method into the R- matrix theory. Within RMHO, it is 
possible to significantly increase the number of states in the daughter nucleus to guarantee 
the convergence of the solution. 

As an example, the RMHO formalism has been applied to the g.s. proton emission from 
141 Ho, in which there have been some experimental hints (e.g., large signature splitting in 
the g.s. rotational band or presence of low-lying 7- vibrational states in the neighboring even- 
even nuclei) for triaxiality. Our calculations show that while the coupling to 7 vibrations can 
in general influence decay characteristics (half-life, branching ratios), in the case of 141 Ho 
the resulting trend is opposite to what has been observed experimentally. From this point of 



view, our results support conclusions drawn in the recent work 



22 1 based on the adiabatic 



particle-rotor approach. An important piece of physics which is still missing in our non- 
adiabatic formalism is the inclusion of quasi-particle pairing. We are currently working on 

n 

incorporating the Hartree-Fock-Bogoliubov couplings |4J] into our model. 
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APPENDIX A: TRIAXIAL COULOMB POTENTIAL 



The Coulomb interaction between the proton and the daughter nuclei is 



r r (r>= / t^W, ( Al ) 
r — r 



where the charge density reads 



P{r) ~ l + ex V [(r-R(n))/a} (A2) 

and the nuclear surface R(Q) is given by Eq. (J2J). To the first order in 02 the Coulomb 
potential is 

Uc(r) = U«(r)| a2=0 + a 2 v4 2) (r). (A3) 
The first term, Vq \ is the Coulomb potential due to an axial symmetric charge density. It 



is given by a simple expression derived in, e.g., Ref. |45j|. The second term, Vq , is of the 
form 

Calculating the derivative of the charge density and taking the limit of the sharp charge 
distribution (a — > 0), one obtains 

V#>(r) = JW^ I MpM22±jkP jy, (A5) 
° V ' V ' J ^/R 2 a (n')+r 2 -2R a (n>)rcos£ 

where 

cos £ = cos 9 cos & + sin 6 sin 9' cos(0 — <t>'). (A6) 
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and R a {Q') = R{Q')\ a 2 =o- This can be reduced to 



vg\v) = R Q p Q C{a Q ,aJ /n 



NT 



f\ dtR 2 Jt)(l - t 2 ) C d<f>' , ^V- 1 (A7) 

where t = cos 8'. The integral over 0' can be calculated analytically, and the final result can 
be expressed in terms of a simple one-dimensional integral: 

V£ 2) (r) = RoPoC(a ,a 2 )J§cos(2<t>) 

<fci£(f)(l - t 2 )^5 {( 4fl2 - &2 )^(^) - 4fl ( fl + & )^(^)} . ( A8 ) 

where 

a = R 2 a (t) + r 2 - 2R a (t)rcos6cos9', (A9) 



b = 2R a (t)r sin 9 sin 9', (A10) 

and 




(All) 



In Eq. ifftBjl . if(/c) and £?( k) are the complete elliptic integral of the first and second kind, 
respectively. 
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